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Abstract 

We review a recent development in theoretical understanding of the quenched 
averaged correlation functions of disordered systems and the logarithmic conformal 
field theory (LCFT) in d-dimensions. The logarithmic conformal field theory is 
the generalization of the conformal field theory when the dilatation operator is not 
diagonal and has the Jordan form. It is discussed that at the random fixed point 
the disordered systems such as random-bond Ising model, Polymer chain, etc. are 
described by LCFT and their correlation functions have logarithmic singularities. 
As an example we discuss in detail the application of LCFT to the problem of 
random-bond Ising model in 2 < d < 4. 
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1 Introduction 



Random systems represent the spatial inhomogeneity where scale invariance is only pre- 
served on average but not for specific disorder realization. The understanding of the role 
played by quenched impurities of the nature of phase transition is one of the significant 
subjects in statistical physics and has attracted a great deal of attention [1]. According 
to the Harris criterion [2] , quenched randomness is a relevant perturbation at the second- 
order critical point for systems of dimension d, when its specific heat exponent a, of the 
pure system is positive. Concerning the effect of randomness on the correlation functions, 
it is known that the presence of randomness induces a logarithmic factor to the correlation 
functions of pure system [3-4] . Theoretical treatment of the quenched disordered systems 
is a non-trivial task in view of the fact that, one has to average the logarithm of the 
partition function over various realization of the disorder in the statistical ensemble and 
therefore find physical quantities [1]. There are two standard methods to perform this 
averaging, the supersymmetry (SUSY) approach, and the well-known replica approach. 
Recently using the replica approach it has been shown by Cardy [5], that the logarithmic 
factor multiplying power law behavior are to be expected in the scaling behavior near 
non-mean field critical points (see also [54]). It is shown also that the results are vahd for 
systems with short-range interactions and in an arbitrary number of dimensions. He con- 
cludes that in the limit of n ^ of replicas the theory possess of a set of fields which are 
degenerate (they have the same scaling dimensions) and finds a pair of fields which form 
a Jordan cell structure for dilatation operator and derives logarithmic operator in such 
disordered systems. It is proved that the quenched disordered theory with Z — 1 can be 
described by logarithmic conformal field theory as well. The logarithmic conformal field 
theories (LCFT) [6-7] are extensions of conventional conformal field theories [8-10], which 
have emerged in recent years in a number of interesting physical problems of WZNW mod- 
els [11-15], supergroups and super-symmetric field theories [16-22] Haldane-Rezzayi state 
in the fractional quantum Hall effect [23-27], multi-fractahty [28], two-dimensional turbu- 
lence [29-31], gravitaitionally dressed theories [32], Polymer and abelian sandpiles [33-35, 
5], String theory and D-brane recoil [36-44], Ads/CFT correspondence [45-52], Seiberg- 
Witten solution to SUSY Yang-Mills theory [53], disordered systems [54-63]. Also the 
material such as Null vectors. Characters, partition functions, fusion rules. Modular In- 
variance, C-theorem, LCFT's with boundary and operator product expansions have been 
discussed in [64-84]. 

The LCFT are characterized by the fact that their dilatation operator Lq arc not diago- 
nalized and admit a Jordan cell structure. The non-trivial mixing between these operators 
leads to logarithmic singularities in their correlation functions. It has been shown [6] that 
the correlator of two fields in such field theories, has a logarithmic singularity as follows, 

< ip[r)ip[r ) >~ |r — r I log \r — r \ + . . . (1) 

In this article we review the conformal field theory (CFT) and logarithmic conformal 
field theory (LCFT), which have appeared in the last decade as a powerful tool for the 
description of the correlation functions of second-order phase transition of pure and dis- 
ordered critical systems near their fixed points. In section 2 and 3 we present a brief 
and self-content review of the CFT and LCFT and their basic tools in d-dimensions. In 
section 4 using the replica method we show that the disordered systems near their fixed 
point can be described by LCFT. As an example we discuss in details the correlation 
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functions of the random-bond Ising model and its connection to LCFT. We give the ex- 
phcit expression of the various types of quenched averaged 2, 3 and 4-point correlation 
functions of the local energy density. We also show that the ratios of these correlation 
functions to the connected ones have specific universal asymptotic and write down these 
universal functions explicitly. 



2 - Conformal Field Theory 

In the following sub-sections wc introduce the necessary techniques and the basic defi- 
nition such as conformal transformation, conformal group, its representation, correlation 
functions. Ward-identities, Virasoro algebra and its representation in 2d conformal field 
theories, etc. 

• Conformal Transformation &; Conformal Group 

Let us start with definition of conformal transformation. 

Definition- 1 : A transformation of coordinates x' — > x is called conformal if it leaves 
the d-dimensional metric g^j^^i, unchanged up to a scalar factor A(x), i.e. 

g'^,{x) = A{x)g^,{x) (2) 

consider an infinitesimal transformation 

xi' ^x'l" ^xi' + ei'ix) (3) 

to order (e): 

where ea;f5 is the covariant derivative of €«. The condition that the transformation be 
conformal gives: 

+ e^;/. = f{x)g^^ (5) 

where 

f{x) = K{x) - 1 

Now we can eliminate f{x). Let us suppose that g^^, — 7)^1, — Diag{l, —1, • • •). Contract 
eq. (5) with r;^^ to get: 

fix) = ^id,e^) = ^id-e) (6) 

This gives us the first relation between f{x) and e. Differentiate eq. (5) w.r.t. and 
after reordering of indices: 

'^df.d^ea = -q^d^f + rj^adu,/ - Vixi^d^f (7) 
Now contract (7) with rj^"^ gives: 

2dhi, = (2 - d)d^f (8) 
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To find final equation for / apply to eq.(8) and to eq.(5) and find: 

(2 - d)d^d,f = ?7^,9V (9) 

now contract above eq. with rj^^ (using rj^^rj^y = d), therefore we find that the / satisfies 
the following equation : 

2(1 - d)d'^f = (10) 

• Conformal Group with d > 2 

For d > 2 eq. (10) reduces to following simple equation, 

a2/ = =^ f{x) ^ A + B^x^ (11) 
Using the relation f{x) = fc^ac": 

6^ — Qjfj^ -\- b^iiX -\- CnuQiX X (12) 

The case 6 = c = ^ infinitesimal translation. Now suppose a = c = 0. 
Substitute = b^i,x^ in eq. (5) gives: 

2 

hf^y + K^, = ^KVt^u (13) 

Therefore b^i, has two part which are antisymmetric and proportional to rim,, i.e. 

b/uu = A77^^ + = -m^^ (14) 

where A is scaling factor. The term which is proportional to 77^,^ is scaling transformation 
and m^y part is an infinitesimal rotation. To understand the meaning of the part c^,^q, x'^x" 
we start with eq.(9) for d > 2: 

(2 - d)d^d,f = r]^,d^f = ^ d^dj = (15) 

or: 

d^dyd -6 = (16) 

dyd-e = -2by ^ e'^ = fe^x ■ X - 2x^'h ■ x (17) 

This is known as special conformal transformation which has combination of inversion 
and translation. For finite transformation: 

x'^ ^x^ + a^ 
x'^ = Xx^ 

x"" = m'^X 
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Generators of conformal group are: (for d > 2) 



D = -ix^di, 

= -i{2x^x'd, - x^d^) (19) 



and the algebra is 



[Kp, Lf,^] = i{r]p^iK^ - Tjp^Ku,) 
[Pp, -f'/ii/] = iivpi^Pv ~ VpvPiJ.) 



[L^u, Lpa-] — i{l]i/pL^a + Tj^crLvp — Tjp^pLy^ — Tji^jL^p) (20) 

Define new generators as: 

J-i,i^= 2^Pi^~'^'^^ 

J-1,0 = D 

Jo,M= l(P, + K,) (21) 

Note that Ja^ = —Jb,a and a, 6 e {— 1, 0, 1, ■ ■ ■ , ci}. Ja^b satisfy the SO{d +1,1) algebra. 
Number of its parameters is \{d-\- 2){d+ 1). Ja,b satisfy: 

[Ja,h, Jc,d] = iiVadJbc + VbcJad " VacJbd " VbdJac) (22) 



• Representation of Conformal Group 

Consider infinitesimal transformation with parameters ou^. We would like to find the 
representation of Tg so that: 



Now define: 
Near the origin: 



usmg: 



one finds, 



^'{x') = (1 - iuj<^Tg)^x) (23) 

V*(0) = S^^m (24) 

exp(ix'*Pc.)L^^ exp(-ix"Pa) 

— Sfiu ~ ^iJ,Pv + ^vP/x (25) 

P^^x) = -id^^x) (26) 



L^M^) = i{x^d,-x,d^)^x) + S^,^x) (27) 
Similarly define the effect of D and on the origin as A and respectively, leads to: 

K^^x) = 

{k^ + 2x^A - x^S^^ - 2ix^x''d^ + tx^d^)<i>{x) 

D^x) = (A - ix''d^)^x) (28) 

Now for finite transformation we can define the quasi-primary field as following. 
Definition 2 : The field $(a;) which under conformal transformation transforms as: 

$(x) ^ $'(x') = Q(x)^^/2$(x) (29) 

is called quasi-primary filed with scaling dimension Aj/2. The relation between fl{x) 
and the Jacobian is, |^| — Qr'^/'^ . 

• Correlation Functions 

Correlation functions of the conformal field theory should transform as following: 

<^\x\)...^\x'^) >= 

0(xi)^i/2 . . . o(xjv)^^/' < $(xi) • • • $(xjv) > (30) 
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The conformal structure implies strong constraints on the correlation functions of the 
theory. These constrains can be found by using the infinitesimal transformation and the 
above transformation. 

One can show that the two, three and four point correlation functions of are given 

by: 

< $i(a;i)$2(a^2) > 



Cl,2 



0, Ai7^A2 ^ ^ 



^liM (S2) 

A1+A2-A3 A2+A3-A1 A1+A3-A2 v^^; 
•^1,2 -^2,3 -^1,3 

/(^^,^i^r-0n^<.4/^-^^-^^- (33) 

xi,3X2,4 a;2,32;i,4 '■' 

where Xij — \xi — Xj\ and A = E^^;^Aj. 

• The Ward Identities 

Consider an infinitesimal transformation as: 

= (1 - ir)''{x)GaMx) (34) 

where GaS are the generators of group and ?7"(a;)'s are infinitesimal functions. Under 
this transformation action changes as: 

SS^-J d^xd^ij^Xi^)) (35) 

where is conserved current corresponding to the transformation (34). On the other 
hand we can find the change of N-point correlation functions under this transformation. 
Defining ^{N) — ^{xi) ■ ■ ■ ^{xn), it can be shown that to order r]{x): 

< 6^N) >=- J d'^xd^ < j^^N) > r]\x). (36) 

Also using explicit expression of transformation (i.e .eq. (34)), : 

AT 

5^N) = -i ^($(xi) • • • Ga^Xi) ■ ■ ■ ^XN))v''{Xi) 
i=l 
N 

= -i / dW{x) E{*(^i) • • • GaHx) ■ ■ ■ ^xn)}5(x - Xi) (37) 

•' i=l 
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Therefore for given (small) rj: 

d 



N 



< j'^{x)^{xi)---^{xN) > 



^tJ25{x- X,) < $(a;i) ■ ■ ■ G„$(a;) ■ ■ ■ $(a;^) > (38) 



i=l 



This is the Ward identity corresponding to current j^{x). 



• Ward identity corresponding to the conformal invariance 

We know that the Stress- Tensor ( T^,y) is conserved current due to the invariance of 5" 
under transformation x'^^ = x^^ + e^^ with constant e's. It's properties arc: 1) T^,^ = Tj,^ 
and 2) d^T^^ = 0. Conserved charges are P" = J d'^'^xT^" . P" as an operator in Hilbert 
space acting as : [Pj,, $] = —id^,^. More generally [Qa, — —iGa^- There is another 
definition of T^^ . Consider the changes in metric as g^i, — > g'^j^ — Qfj^v + Sg^i,, under this 
transformation action S transforms as : 

5S^ j d'x^T^^Sg,, (39) 

This enable us to find more restrictions of T^'' for conformal transformation. Suppose 
that theory possess Wyle symmetry so that: 

g^^{x) A{x)g^^{x) (40) 

For infinitesimal transformation, g^j^^, — > g^^, + u!{x)gi^,^, or Sg/^i, — uj{x)gnv. Substitute 
this result in eq.(39) : 



= J d'^x^T^''uj{x)g^ 
J d^'x^THuix) = 



^ 7^ = (41) 
This means that the stress-tensor is tracelcss. 

Now we can write the Ward identities (WI) due to conformal invariance. For invariance 
under translation: 

N Q 

d, < Tli^N) >= - ^i)^ < > (42) 

i=l '^■^i 

for Lorantz invariance ( its current is j^^"' — T^'^x°^ — T^x'^ and its generators is 
Lfxv — Sfxi, -\- i{x^di, — Xj^dn)): 

< {Ti'^'x'^ - Ti"^x'')^{N) >= 
- f: 5(x - x,){x\dt - - iSr) < ^N) > (43) 



1=1 
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Using the first identity: 



N 

iY,S{x- xO-^r < HN) > . (44) 



WI corresponds to the scahng invariance 



y: s{x - + A,) < m) >} (45) 

i=l ^^i 

Also using the first identity: 

TV 

< T;^$(Ar) >= - < *(^) > ■ (46) 

i=l 

• Conformal invariance in 2- dimensions 
Suppose X = {x^,x^) and gf^i, — Diag{l, 1). Using eq.(5) we have: 

di^e^ + d^e^ = '^9nu{d ■ e) (47) 
For V ^ and = = 1 we have the following equations 

d\ei = -92ei 

9iei = 9262 (48) 
Define z — x^ ^ ix^, z — x'^ — ix^, and (Ps — dzdz so that : 

9zz = Qzz^O 

Qz-z = gzz = ^ (49) 

now for w = ei + ie2, the eqs. (48), reduce to the condition dzw{z,z) = 0. This 
shows that the group of conformal transformations in two dimensions is isomorphic to the 
(infinite-dimensional) group of arbitrary analytic coordinate transformation z —>■ w{z) 
and z — >• iv{z). 

The Mobius transformation is a subset of holomorphic transformation which has the 
following expression: 

^^/(^) = ^i±^ ad-bc^l. (50) 

^ ^ cz + d ^ ^ 

where a, b, c and d are complex numbers. One can show that the two transforma- 
tions fi and /2 gives, /1/2 = /, so that the parameters of fi and /2 be Oi, 61, Ci, di and 
02, b2, C2, d2, respectively, the parameters of /a will be 03 = 0102+6102, 63 = O162+&1CI2, C3 — 
01O2 + 01C2, d3 — C162 + ai02. 
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• Holomorphic form of conformal ward identity 



In 2D we use the two-dimensional expression for generator of spin of i-th field as Si^^ = 
Sie^y where e^^, = —e^^. Therefore conformal Ward Identities reduces to the following 
form: 

d 



N f) 

i=l ^-^j 

N 

1=1 

N 

< r^^(x)$(7V) >= - E ^(^ - < HN) > (51) 



1=1 

Now using the following identity: 



S'ix) = -d,- = -dA (52) 
n Z TT z 



we can write the conformal WI as: 

N ^ 

27rd, < T,,^N) > +27rd, < T.-^N) > 

N ^ 

i=i ^ 

2 < r,,$(7V) > +2 < T,MN) > 

N 

= -5^5(x-xOA, <$(7V) >, 

i=l 

-2 < T,-MN) > +2 < T,MN) > 

N 



-^5{x-Xi)si<^N) >, (53) 



1=1 



Using the last equations: 

27r < T,MN) > 



N ^ 



i=i ^-W' 
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27r < T,,$(iV) > 

TV ^ 

^ - Y.9z{^^)hi < $(iV) >, (54) 



1=1 



where /i = i(A + s) and h = |(A — s). Using the above equation we can rewrite the 
eq.(53) as: 

d, < f $(iV) >= 

MY. <m)> + < m) >} (55) 

where T = -27rT^, and f = -27rT^^. Therefore: 



< T(z)$(Ar) 

g{_J_a„,<*(^)>+^<*W >, + ... (56) 

Also same equation holds for T (replace z with z). This gives us the Operator Product 
of T and For example consider = 1 and find that : 

T{z)^w) = -^d^^w) + . ^ + ■ ■ • (57) 

z — w [z — wy 

In the next section we will start from the OPE of T and $ and introduce the LCFT. 
• Correlation Functions in Two-Dimensions 



Let us now find S^^{x). Using the eq.(36): 

6,<^N)>= I (fxd^ < T>"'{x)e^{x)^N) > 
Jm 

I- (f{-dz < T''e,^N) > +dz < T''e,^N) >} (58) 
2 Jc 



or 



6,, < $(iV) >= — — (fdz< T{z)<^{N) > +C.C (59) 
27ri Jc 
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Using the operator product of T{z) and ^{w): 



N 



< $(iV) >= - ^{e{wi)d^^ + d^^e{wi)hi) < $(iV) > (60) 

i=l 

Therefore for holomorphic part: 

5,^z) = -ed.^z) - h^z)d,e (61) 
For infinitesimal transformation e{z) — a + bz + cz^: 

N 

1=1 

N 

Y^iwid^, + hi) < $(iV) >= 

i=l 

N 

^{w^d^, + 2wihi) < $(iV) >= (62) 
1=1 

One can solve the above equation and find that: 

< ^i{zi, zi)^2iz2, Z2) >= 5h,M 7 (63) 

[Zi - Z2Y'^{Zi - Z2Y'^ 



< $i(xi)*2(a;2)*3(2;3) 



(64) 



hi+h2—h3 h2+h3—hi hi+hz—h^ 
•^1,2 -^2,3 -^1,3 



< ^l{Xi)^2{x2)^z{.X3)^i{Xi) > = 

/(^^,^i^,---)nf<,4f-'^-'^^- (65) 

2^1,33^2,4 2^2,32:1,4 

where h — Z^^^i hi. Note that in four point function unknown function / has only one 
crossing ratio, because, if one define n — ^i^^^ therefore will find that, :5i^2£M — J_ 

oil I Z2, 3^1,4 ' ' 2:1,422,3 1— »7 

= 1 - 77. So / is function of rj and fj. 

In the end of this sub-section let us introduce the central charge c via the OPE of T 
and T. For instance consider the following actions in two dimensions; 

Si^\g j d^xd^(f>d^<p 

S2^\g j d^xi^+^^^d^^ (66) 
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it can be shown that : 

rwrM = ^.JIH (67) 

[z — w)^ [z — wy z — w 

where c—1 and c— 1/2 for 5"! and S2, respectively. This is just the definition of central 
charge c. 

• Transformation of Stress- Tensor 

According to definition of stress-tensor: 

1 



and for T{w): 



5,^ = -— i e{z)dzT{z)^ (68) 
S,T{w) = --^ i e{z)dzT{z)T{w) 

= -^cdle{w) - 2T{w)d^e{w) - e{w)d^T{w) (69) 

This gives the transformation of T under infinitesimal transformation z ^ z-\- e{z). For 
finite transformation { z ^W{z)): 

T'{W) = (^^)-'[T{z)-^{W-z}] (70) 

where {W; z} = '^^^^.^ 3/2( ^^|^^^^ )^, which is called the Schwartzian derivative. 

As an example one can check that for Mobius transformation W = (with ad — bc= 1) 
we have {W; z} = 0. 

Also the map W{z) — L/27r\nz maps the plane to the strip with periodic boundary 
conditions. One can show that 1) < Tpiane >— and 2) < Tcyi >— 

• The Virasoro Algebra & 
Its Representation 

Expanding the stress-tensor in terms of a Laurent series as: 

r(^)=X:^"""'^n (71) 



with 



K = — i dzz'^+^Tiz) (72) 



27ri / 

Using the above expansion and operator product expansion (OPE) of TT: 



dw r dz 



2m 
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c/2 ^ 2T{w) ^ dT{w)^ 



(z — w)^ (z — wY z — w 
J 27ri 12 

+2(n + l)w"+™+^T(w) + w''+"'+^dT{w)) 

[Ln, Lm] = Y2^i'^^ - l)<^n+m,0 + {u - m)Ln+m (73) 

These commutation relations are known as Virasoro Algebra. We intend to create a 
representation of this algebra in terms of states, so let |0 > be a vacuum state in our 
field theory such that L„|0 >= for n = 0,1,-1. Define the state \h >= $(0)|0 > 
for any primary field $ of conformal weight h. It can be shown that Lo\h >— h\h >, 
so that \h > are eigenvectors of Lq. Also we have Ln\h >— for n > 0. We call 
these states highest — weight states in the representation. We can obtain all states in 
the representation by applying a sequence of L_„ operators to a highest weight state. 
Alternatively we can define fields: 

= L^k, ■ ■■L-kMz) (74) 

with 

^ ^ . . r dz T{z)^{w) 



where /ci H \-kn = N. The fields $-'=i'-'-'="(^) are known as the descendent of field 

^{z) at level N. 

For example one can check that the stress-tensor is the descendent of identity operator 7, 
i.e. L_2/ = T(0). 

Now define null vector |(^ > as being a linear combination of states of the same level. 
Similar to primary state it is it must be stable under the operation of L„ and Lq. 

Ln\i > =0, n > 



Lo|e> ={h + N)\^> (76) 

For instance we can show that the state (L_2 + aL^j)$ with a — — 2{2h+i) ^ 
state in the level N = 2. The existence of null vectors gives us additional differential 
equation for determining of unknown function f[rj) in four-point correlation functions 
which contains the field $. 
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3 -Logarithmic Conformal Field Theory 

In an ordinary conformal field theory primary fields are the highest weights of the 
representations of the Virasoro algebra. The operator product expansion that defines a 
primary field ^{w,w) is 

T{z)^i{w,w) = ,^ Mw:W) + , ^ d^^,{w,w) (77) 

[z — wy [z — w) 

T{z)^i{w, w) = w) + , ^ _. d^^i{w, w) (78) 

\Z — W) \Z — W) 



where T{z) := T^z^z) and T{z) := T^z^z). 

The primary fields are those which transform under z — > f{z) and 
z f{z) as: 

^i{z,z) '^i{z,zy = 

i%^)'''i%^f'Hf-\z),f-\-z)) (79) 

It can be shown that the OPE of T and $ is equivalent to the commutation relation 
between L^s and $. Let us introduce the radial quantization which is defined as: 

R(^i(z)^2iw)) = <^ ^; \J", ( (80) 

Now suppose we have two function a(z) and h{z) (are holomorphic) and consider the inte- 
gral §^ dza{z)b{z) . To have a well defined expression for this integral we should consider 
the radial ordered of fields. We consider two contours Ci and C2 so that the ci is the circle 
path ( anti-clockwise) with center in the origin and its radius is r = |w| + e and C2 is 
similar path as Ci ( but clockwise) with r — \w\ — e. Therefore, 

dza{z)b{z) — i dza{z)b{w) — i dzb{w)a{z) (81) 

Jc J C\ J C2 

Define A = §dza{z), therefore, 

£dza{z)b{z) ^ [A,b{w)] (82) 

Prom eq. (72) one finds, 

[Ln, Hw)] = ^ / dzz^+'T{z)^w) (83) 

Using the OPE of T and $ 

[L„, ^i{z)] = z'^+'dz^i + (n + 1);2"A,$, (84) 
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One can regard Aj's as the diagonal elements of a diagonal matrix D, 

[L„, ^i{z)] = z^^'d.^i + {n + l)z''Dl^j (85) 

One can however, extend the above relation for any matrix D, which is not necessarily 
diagonal. This new representation of L„ also satisfies the Virasoro algebra for any arbi- 
trary matrix D. Because we have not altered the first term in the right hand side of the 
equation (85), this is still a conformal transformation. 

By a suitable change of basis, one can make D diagonal or Jordanian. If it becomes 
diagonal, the field theory is nothing but the ordinary conformal field theory. The general 
case is that there are some Jordanian blocks in the matrix D. The latter is the case of a 
LCFT [7]. 

Here, there arise some other fields which do not transform like ordinary primary fields, 
and are called logarithmic operators. For the simplest case, consider a two-dimensional 
Jordan cell. 

The fields $ and satisfy 

[L„, $(^)] = z"+^(9,$ + (n + 1)^"A$ (86) 

and 

^(z)] = + (n + l)z''A^ + {n+ 1)/'$, (87) 

and they transform as below 

m - i%^)^'^u-\^)) (88) 

m - (^)^[*(/-^(^)) + ^og{^^^Mf-\z)] (89) 

Note that we have considered only the chiral fields. The logarithmic fields, however 
cannot be factorized to the left- and right-handed fields. For simplicity we derive the 
results for chiral fields. The corresponding results for full fields are simply obtained by 
changing 

^ z'^z^ (90) 

and 

log 2; ^ log 1^1^ (91) 

Now compare at the relations (86,87) and (88,89); one can assume the field ^ as the 
derivation of the field $ with respect to its conformal weight, A. Now let us consider the 
action of Mobius generators (Lq, L±) on the correlation functions. Whenever the field \E' 
is absent, the form of the correlators is the same as ordinary conformal field theory. By 
the term form we mean that some of the constants which cannot be determined in the 
ordinary conformal field theory may be fixed in the latter case. Now we want to compute 
correlators containing the field \E'. At first we should compute the two-point functions. 
The two-point functions of the field <l> is as below 

< ^zMw) >= ^— ^ (92) 
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In the ordinary conformal field theory the constant c cannot be determined only with 
assuming conformal invariance; to obtain it, one should know for example the stress- 
energy tensor, although for c 7^ one can set it equal to one by renormalizing the field. 

Assuming the conformal invariance of the two-point function < \E'(z) $(■«;) >, means 
that acting the set {Lo,L±i} on the correlator yields zero. Action of L_i ensures that 
the correlator depends only on the z — w. the relations for L+i and Lq are as below 

[z'^d^ + w^dyj + 2A(z + w)] < > + 

2z < ^{z)^{w) >= (93) 



[zd^ + wd^ + 2 A] < > + 

< ^{z)^{w) >= (94) 

Consistency of these two equations for any z and w, fixes c to be zero. Then, solving the 
above equation for < ^(^) $(■«;) > leads to 

< ^{z)^{w) >= 0, 

< ^zMw) >= j^^^ (95) 

Now assuming the conformal invariance of the two-point function < ^(^)^('u;) >, gives 
us a set of partial differential equation. Solving them, we obtain 

< ^{z)^{w) >^ _^ [b - 2a \og{z - w)] (96) 

Now we extend the above results to the case where Jordanian block is n -|- 1-dimensional. 
So there is n -|- 1 fields with the same weight A. 

[L„, ^i{z)] = + {n + l)z"A$i + {n + l)^"$i-i, (97) 

where $_i = 0. 

All we use is the conformal invariance of the theory. Prom the above fields, only $0 is 
primary. Acting L_i on any two-point function of these fields, shows that 

<<^>i{z)<^j{w)>= fij{z-w). (98) 

Acting Lq and L+i, leads to 

< [Lo, $i(^)$,(0)] >= (zd, + 2A) < M^)^A^) > 

+ < $,_i(^)$,(0) > + < $,(z)$,_i(0) >= (99) 

< [L+i,$,(z)$,(0)] >= {z% + 2zA) < $i(^)$,(0) > 

+2z < $i_i(^)$j(0) >= 0. (100) 
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Then it is easy to see that 

< $,_i(^)$,(0) >=< $,(^)$,_i(0) > . (101) 
Using $_i = and the above equation, gives us the following two-point functions. 

< ^i{z)^j{w) >= for i+j<n (102) 

Now solving the Ward identities for < ^o{z)^n{w) > among with the relation (101), leads 
to 

< ^i{z)^n-iH >=< %{z)^n{w) >= ao{z - H"'^- (103) 
The form of the correlation function < $i(2;)$„(u') > is as below 

< $i(^)$„(w) >= {z - w)-2^[ai + 6i log(^ - w)], (104) 
but the conformal invariance fixes bi to be equal to — 2ao. So 

< ^i{z)^n+l-i{'w) > = < ^i{z)^n{'w) > = 

{z - w)-2^[ai - 2ao log(;2 - w)] for i > (105) 

Repeating this procedure for the two-point functions of the other fields $i with 
and knowing that they are in the following form 

< M^)^nM >= (z-wy^^j^^'iA^^siz-w))^, (106) 

3=0 



gives 



or 



So 



^ia,,(log(^ - w)y-' + 2 a^-l,J{log{z - w)y = (107) 

j=l j=0 



or 



(-2)^'+' _ (-2)^'+! 

j——^ai-j-i,o -: j-—^ai-j-i (108) 

< Mz)K{w) >= {z - w)-2^ ^-^ai_,{log{z - w)y, (109) 

j=o J- 

and also we have 

< <^>i{z)<^>k{w) >= 

< ^i+k-n{z)^n{w) > for i + k>n. (110) 
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So for the case of n logarithmic field, we found all the two point functions. The interesting 
points are 

i) some of the two-point functions become zero. 

ii) some of the two-point functions are logarithmic, and the highest power of the logarithm, 
which occurs in the < >, is n. 

The most general case is the case where there is more than one Jordanian block in the 
matrix D ,or in other words, there is more than one set of logarithmic operators. The 
dimension of these blocks may be equal or not equal. Using the same procedure, one can 
find that 



^((^-^)-^^EItfc^'^^«n-.[log(^-^)]^ ^ + J>n (111) 
i 0, i + j < n 

where / and J label the Jordan cells, n — max{n/, rij} and rii and rij are the dimensions 
of the corresponding Jordan cells. 

Also note that the conformal dimensions of the cells I and J must be equal, otherwise 
the two-point functions are trivially zero. 

Now we want to consider the three-point functions of logarithmic fields. The simplest 
case is the case where, besides $, only one extra logarithmic field \E' exists in the theory. 
The three-point functions of the fields $ are the same as ordinary conformal field theory. 

= 777T^=^«/(^i'^2,e3), (112) 
where ^ 

If one acts the set {Lo,-L±i} on the three-point function < ^ {zi)^{z2)^{z3) >■— 
B(zi, Z2, Z3), the result is an inhomogeneous partial differential equation for B{zi, Z2, Z3) 
where the inhomogeneous part is 74(^1, ^2, ^3)- So the form of B{zi, Z2, Z3) should be as 
below, 

Biz^,Z2,Zs) = [b + Y.b^^Og^^]f {^1,^2,^3). (US) 

Solving the above mentioned differential equations, we find the parameters bi to be 

bi = -62 = -bs = a. (114) 

The final result is 

< ^{z^Mz2Mzs) >= [6 + alog-^]/(ei,6,e3) (115) 

?2?3 

If there are two fields ^ in the three-point function, one can write it in the following form 

< *(^i)*(^2)*(^3) > = 
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[C+ ^QlogCi + Ec^.l0ge^l0gei]/(el,6,e3)• (116) 
ij 

Again the Ward identities can be used to determine the parameters q and q^, 

< *(^i)*(z2)$(^3) >= 

[c - 261026 + a[{-^r + (log 6)1/(6, 6, 6)- (117) 
Finally, for the correlator of three ^'s we use 

< *(Z1)*(Z2)*(^3) >= [d + diDi + ^2^2+ 

+ ^3^3 + ^/ ^^^^ + 4^?]/(ei, 6, 6) (118) 

where 

A := log(666) (119) 

D2 := log 6 log 6 + log 6 log 6 + log 6 log 6 (120) 

A3 = log 6 log 6 log 6- (121) 

This is the most general symmetric up to third power logarithmic function of the relative 
positions. Using the Ward identities, this three-point function is calculated to be 

< ^'(zi)*(z2)*(z3) >= [d - cDi + AbD2 - hDl 

+8aDs - AaD,D2 + aDl]f{^,, 6, 6) (122) 

Now there is a simple way to obtain these correlators. Remember of the relation between 
^{z) and ^(z) 

^(z) = ^^z). (123) 

Consider any three-point function which contains the field ^P. This correlator is related 
to another correlator which has a $ instead od ^ according to 

< ^(^1)^(^2)5(^3) > 

= ^ < $(^1)^(^2)5(^3) >, (124) 

To be more exact, the left hand side satisfies the Ward identities if the right hand side 
does so. 

But the three-point function for ordinary fields are known. So it is enough to differentiate 
it with respect to the weight A. Obviously, a logarithmic term appears in the result. 
In this way one can easily obtain the above three-point functions. In fact instead of 
solving certain partial differential equations, one can easily differentiate with respect to 
the conformal weight. This method can also be used when there are n logarithmic fields. 
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To obtain the three- point function containing the field $i, one should write the three- 
point function, which contains the field $0, and then differentiate it i times with respect 
to A. 

Note that in the first three-point function, there may be more than one field with the 
same conformal weight A. Then one must treat the conformal weights to be independent 
variables, differentiate with respect to one of them, and finally put them equal to their 
appropriate value. Second, there are some constants, or unknown functions in the case 
of more than three-point functions, in any correlator. In differentiation with respect to a 
conformal weight, one must treat these formally as functions of the conformal weight as 
well. 

As an example consider 



< $(zi)$(z2)^(-Z3) > = 

a 



(125) 



Differentiate with respect to Ai, and then put Ai = A2 = A3, and = b. This is 
(115). This method can be used for any n-point function: 



< <^>i{Zi)---A{Zn_i)B{Zn) > = 

^ < ^z,) ■ ■ ■ A{zn-i)B{z^) >, (126) 

provided one treats the constants and functions of the correlator as functions of the 
conformal weight. 

• Logarithmic Conformal Field Theory With Continuous Weights 

In the previous subsection assuming conformal invariance we have explicitly calculated 
two- and three-point functions for the case of more than one set of logarithmic fields when 
their conformal weights belong to a discrete set. Regarding logarithmic fields formally as 
derivations of ordinary fields with respect to their conformal dimension, we have calculated 
n-point functions containing logarithmic fields in terms of those of ordinary fields. Here, 
we want to consider logarithmic conformal field theories with continuous weights [72]. 
The simplest example of such theories is the free field theory. 

In the last part it is shown that if there are quasi-primary fields in a conformal field 
theory, it causes logarithmic terms in the correlators of the theory. By quasi-primary 
fields, it is meant a family of operators satisfying 



[L„, $(^)(^)] = z''+^d,^^'\z) + {n + 1)^"A$(^)(^) + (n + l)z''A<l>^'-^\z), (127) 

where A is the conformal weight of the family. Among the fields the field 

is primary. It was shown that one can interpret the fields formally, as the j-th 

derivative of a field with respect to the conformal weight: 

^^^'H^) = -T^<^^°H^), (128) 
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and use this to calculate the correlators containing in terms of those containing 
$(0) only. The transformation relation (127), and the symmetry of the theory under 
the transformations generated by L±i and Lq, were also exploited to obtain two-point 
functions for the case where conformal weights belong to a discrete set. There were two 
features in two point functions. First, for two families $i and $2; consisting of ni + 1 
and 712 + 1 members, respectively, it was shown that the correlator < ^[^^^2 ^ > is zero 
unless i + j > max(ni,n2). (It is understood that the conformal weights of these two 
families are equal. Otherwise, the above correlators are zero.) Another point was that 
one could not use the derivation process with respect to the conformal weights to obtain 
the two-point functions of these families from < >, since the correlators contain a 

multiplicative term 5ai,A2) which can not be differentiated with respect to the conformal 
weight. 

Now, suppose that the set of conformal weights of the theory is a continuous subset of 
the real numbers. First, reconsider the arguments resulted to the fact that < ^^^^2^ > 
is equal to zero for i + j > max(ni,n2). These came from the symmetry of the theory 
under the action of L±i and Lq. Symmetry under the action of L_i results in 

< ^f\z)<^^^\w) >=< ^f\z - w)^^^\0) >=: A'^{z - w). (129) 

We also have 

< [Lo, $?(^)*^^'^(0)] {zd + Ai + A2) A'^{z) + A'-''^{z) + A''^-\z) = 0, (130) 

and 

< [Li, ^f\z)^^^\0)] >= {z^d + 2zAi)A'^{z) + 2zA'-^'^{z) = 0. (131) 
These show that 

(Ai - A2) A'^{z) + A'-^^^{z) - A'^^-\z) = 0. (132) 

If Ai 7^ A2, it is easily seen, through a recursive calculation, that A^^^s are all equal to 
zero. This shows that the support of these correlators, as distribution of Ai and A2, is 
Ai — A2 = 0. So, one can use the ansatz, 

AHz) = T.^ki^)S^'\^i-^2). (133) 

k>0 

Inserting this in (132), and using x5^''~^^\x) — —{k + l)5^''\x), it is seen that 

E[-(^ + iXi(^) + ^k-''{z) - A^'-\z)]5^'\A^ - A2) = 0, (134) 

fc>0 

or 

(k + l)Ati(^) = A'-^'^z) - A'^-'iz), k > (135) 
This equation is readily solved: 

1=0 

where Aq^s remain arbitrary. Also note that A^^^s with a negative index are zero. We 
now put (133) in (130). This gives 

{zd + Ai + A2)^;^(;^) + Al-'^z) + Ai^-\z) = 0, (137) 
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Using (136), it is readily seen that it is sufficient to write (137) only for /c = . This gives 
{zd + Ai + A2)Ai^{z) + Ai-''\z) + Ai^^-\z) = 0. (138) 

Putting the ansatz 

^^^■(^) = ^-(^^+^^)X:«^(ln^)'" (139) 

m=0 

in (138), one arrives at 

(m + l)aZ+i + cC,''' + c^-' = 0, (140) 

the solution to which is 

/ 1 \ m m 

Prom this 



«™ = EiTWo-""^''"- (141) 



^i{z) = ;.-(^^+^^) y: i^^^r^-;^^ E(^)«o-'"^^''"^ (142) 

and 



m=0 s=0 



1 k i+j-k ni 

^ii^) = bE(-l)'(f) E (ln^)™^E(r)«o-'-"^^'+^"-'-1^-^^^+^^^- (143) 



So we have 



A^^(.) = .-(^^+^^)E<^^'^(^i-^2)[^E(-iy(f) 



fe>o 1=0 



i+'j-k ( -j^-jm m 

m 



(ln^)™^^^(-)a^-*^— +'+^'^-'-^], (144) 



m=0 s=0 

or 



^*i(^) = ^-(^1+^2) y L4f|I^a'-™-«-*(ln^)''+^5(P+'^)(Ai - A2), (145) 
where 

q;^^' := ot^. (146) 

These constants, defined for nonnegative values of i and j, are arbitrary and not deter- 
mined from the conformal invariance only. 

Now differentiate (145) formally with respect to Ai. In this process, a'-^'s are also 
assumed to be functions of Ai and A2. This leads to 



dA'^{z) 

— z 



( ■\\q+r+s plrui—p—T,j—q—s 

"^^^^^^^ E ^-i^y:r{-^K (in.)^'+^-5(^+^nAi-A2) 



or 



9Ai p\q\T\s\ ' OA, 

+ Q;^-™-«-^[(lnz)^+"5^+«+^)(Ai - A2) - (lnz)^+"+M(f+«)(Ai - A2)I>,47) 



^^^'(^) _ ,-(A.+A.) J2 ^ zY+^5^^+'^\ A^-A2) 



dAi p^^^ p\q\r\s\ 
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[{p + rja^-f-'-'^-"-* + — ]. (148) 



Comparing this with A^^^'\ it is easily seen that 



provided 



f)ryi-p-r,3-<l-s 

^ + 1 _ p _ r)Q;^+^-™-^-^ (150) 

Note, however, that the left hand side of (150) is just a formal differentiation. That is, 
the functional dependence of a'-^ 's on Ai and is not known, and their derivative is just 
another constant. Repeating this procedure for A2, we finally arrive at 

and 

These relations mean that one can start from which is simply 

A^\z) - ^-(^i+^^)5(Ai - A2)a°°, (153) 

and differentiate it with respect to Ai and A2, to obtain A^K In each differentiation, some 
new constants appear, which are denoted by a'^'s but with higher indices. Note also that 
the definition is self-consistent. So that this formal differentiation process is well-defined. 

One can use this two-point functions to calculate the one-point functions of the theory. 
We simply put $5°^ = 1. So, A2 = 0, 

< $(o)(^) >=/?05(A), (154) 

and 

i an—k 

<$«(^)>=E V'^'(A), (155) 
fe=o 

where 

(156) 

The more than two-point function are calculated exactly the same as in previous subsec- 
tion. 

• The Coulomb-gas model as an example of LCFT's 

As an explicit example of the general formulation of the previous section, consider the 
Coulomb-gas model characterized by the action [8-9] 

-5 = ^ / d''x^g\-g^\d^^){d,^) + iQmi (157) 
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where $ is a real scalar field, Q is the charge of the theory, R is the scalar curvature of 
the surface and the surface itself is of a spherical topology, and is everywhere fiat except 
at a single point. 

Defining the stress tensor as 



:= -i^i^. (158) 



it is readily seen that 



and 

T{z) := T,,(z) = -(90)2 - iQd^, (160) 
where in the last relation the equation of motion has been used to write 

$(^,z) =0(^) + 0(^). (161) 

It is well known that this theory is conformal, with the central charge 

c = 1 - 6Q2 (162) 

There are, however, some features which need more care in our later calculations. First, 
this theory can not be normalized so that the expectation value of the unit operator 
become unity. In fact, using as the integration measure, it is seen that 

< 1 >oc 5(Q) (163) 
one can, at most, normalize this so that 

< 1 >= 5(Q) (164) 

Second, has a 2;-independent part, which we denote it by 0o- The expectation value of 
00 is not zero. In fact, from the action (157), 

< >=< 00 y 90000 exp(2ig0o), (165) 

where N is determined from (164) and 

< 1 >= ^ / 900 exp(2^g0o). (166) 
This shows that A'"(0) = tt, and 



More generally 
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Third, the normal ordering procedure is defined as following. One can write 

(/)(^) = 0o + 0+(-2)+0-(^), (169) 

where < O|0_(^) = 0, 0+(^)|O >= 0, and 

[0o,0±] = O. (170) 

The normal ordering is so that one puts all '-' parts at the left of all '+' parts. It is then 
seen that 

<: m :>=< /(0o) > (171) 

Here, the dependence of / on in the left hand side may be quite complicated; even / 
can depend on the values of 4> at different points. In the right hand side, however, one 
simply changes 0(^) —>■ 4>o- 

Now consider the two point function. Prom the equation of motion, we have 

< (t){z)(t){w) >= -]- \n{z - w) < 1 > +6 (172) 

we also have 

<: 0(^)0(^) :>=< 4>l >= -^^[^'^(Q)] (173) 

Note that there is an arbitrary term in (172), due to the ultraviolet divergence of the 
theory. One can use this arbitrariness, combined with the arbitrariness in N[Q), to 
redefine the theory as 

0(^)0(«;) =: Hz - w)+ : <l>{z)<l>{w) :, (174) 

and 

these relations, combined with (171) arc sufficient to obtain all of the correlators. One 
can, in addition, use (160) (in normal ordered form) to arrive at 

T(.)0(H = ^-^i^ + r.t., (176) 

z — w [z — wy 

and 

d^T 2T{w) , (l-6Q^)/2 

T{z)T{w) = + — — . 177 

z — w [z — wy yz — wy 

Eq. (176) can be written in the form 

[L,, <P{z)] = z'^+'d<P - Y(n + l)z\ (178) 

This shows that the operators (f) and 1 are a pair of logarithmic operators with A = (in 
the sense of first part of section 3). Note that < 1 > is equal to 5{Q). One can easily 
show that 

r(.) : : = ^ ^ ■ - + : : +.t., (179) 

z — w [z — wy 
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which shows that : e'""^ : is a primary field with 

a{a + 2Q) 



(180) 



To this field, however, there corresponds a quasi conformal family (pre-logarithmic oper- 
ators [55]), whose members are obtained by explicit differentiation with respect to a (a 
is not the conformal weight but is a function of it): 

W(") : = : (/.V""^ : = (-i)"7^ : e^""^ : . (181) 

To calculate the correlators of VF's, it is sufficient to calculate < • ■ ■ W'^^ >. 
One has, using Wick's theorem and (174), 

n^_i : e^"'^'^^^^^ : = e^^^^i^'<J^'="'"^'^°^^*~^^'^ : e*^J^=i "^'^^^^^ : (182) 
Prom this using (171) and (174), we have 

1 ^ 

= [^i<i<j<k{zi-zj)^]S{Q + -Y,aj). (183) 

Obviously, differentiating with respect to any Oj, leads to logarithmic terms for the 
correlators consisting of logarithmic fields W^^\ The power of logarithmic terms is equal 
to the sum of superscripts of the fields VF^"^. 

• Logarithmic Conformal Field Theory in d-dimensions 

In the previous subsection we discussed the LCFT in 2-dimensions. Generahzation for 
arbitrary dimension d has been given in [68] . Wc have dealt with two dimensional confor- 
mal field theory relying heavily on the underlying Virasoro algebra, and have described 
how the appearance of logarithmic singularities is related to the modification of the rep- 
resentation of the Virasoro algebra. In this subsection we try to understand LCFT's in 
the context of d-dimensional conformal invariance. 

As is well known, one of the basic assumptions of conformal field theory is the existence of 
a family of operators, called scaling fields, which transform under scaling S : x ^ x' — Xx, 
simply as follows: 

(j){x) (p'ix') = \-'^(j)ix), (184) 

where A is the scaling weight of 4>{x). It is also assumed that under the conformal group, 
such fields transform as, 

0(x)^0V)=ll^ll^0(^), (185) 
(px 

where d is the dimension of space and || ^ || is the Jacobian of the transformation. Equa- 
tion (185) which encompasses eq. (184) defines the transformation of the quasi-primary 
fields. For future use we note that the Jacobian equals A'' for scaling transformation and 
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X 



I ^'^ for the Inversion transformation / : x 



X 



being unity for the other 



elements of the conformal group. Combination of (185) with the definition of symmetry 
of the correlation functions, i.e.: 



< ■ ■ ■ (I>'n{^'n) >=< M^'i) ■ ■ ■ 0Jv(a:^) >, 



(186) 



allows one to determine the two and the three point functions up to a constant and the 
four point function up to a function of the cross ratio. 

It's precisely the assumption that scaling fields constitute irreducible representations of the 
scaling transformation, which imposes power law singularity on the correlation functions. 
As we will see, relaxing this assumption, one naturally arrives at logarithmic singularities. 
It also leads to many other peculiarities, in the relation between correlation functions. To 
begin with, we consider a multiplet of fields. 



^2 



(187) 



V 0n / 

and note that under scaling a; — > Ax, the most general form of the transformation of $ is, 

$(x) ^ $'(x') = X^'^x) (188) 
where T' is an arbitrary matrix. More generally, we replace (188) by, 

<^{x) ^ <^'{x') =11 ^ f $(x). (189) 

where T is an n x n arbitrary matrix. When T is diagonalizable, one arrives at ordinary 
scaling fields by redefining $, so that all the fields transform as 1-dimensional represen- 
tation. Otherwise, following [68] we assume that T has Jordan form. 



' d 
1 



V 



••• 


\ 


-A 





d 


1 ■•. 





... 1 


d 1 



(190) 



Rewriting T as -^1 + J, where Jjj = ^j-ij, eq. (189) can be written in the form. 



d>x 



$(x) ^ $'(a;') =11 ^ 11^ ^x^ix), 
where A^. =|| ^ ||"^ is a lower triangular matrix of the form, 

{In II f 



[191] 



i. e. for iV = 2 we have, 
<P[{x') 

02 (x) 



= 1, 



(192) 



ox 



-A 

d 



(px 



(f)x' 

d>x 



(In I 
28 



(j)x' 
(hx 



(pi{x) + 02 (a;)). 



(193) 



An important point is that the top field (j)i{x) always transform as an ordinary quasi- 
primary field. A most curious property of the transformation (191) is that each field (pk+i 
transforms as if it is a formal derivative of (f)k with respect to 

This formal relation which determines the transformation of all the fields of a Jordan cell 
from that of the top field 0i, essentially means that with due care, one can determine the 
correlation functions of the lower fields from those of the ordinary top fields simply by 
formal differentiation. Also one can find the two-point, three-point correlation functions 
of fields for Jordan -cell of rank r. For example we already know by standard arguments 
that the two point function of the top fields (pa and belonging to two different Jordan 
cells {Aa,n) and (A^,m) vanishes, i.e.: 

\\ X — y 

Due to the observation (194), it follows that the two point function of all the fields of two 
different Jordan cells with respect to each other vanish. Therefore we can calculate the 
two point function of the fields within the same Jordan cell. As we will see logarithmic 
conformal symmetry gives many interesting and unexpected results in this case. Let's 
denote the matrix of two point functions < (f)i{x)(()j{y) > for all 0j, (pj e (A, n) by G{x, y), 
then from rotation and translation symmetries, this matrix should depends only on || 
X — y II . From scaling symmetry and using (185) and (186), we have, 

AG'(||x-y||)A* = A2^G'(A||x-y||), (196) 

where A = A*^"^, and from inversion symmetry, we have, 

A.G(|| X - y ||)A; =11 X - y \\-'^" G(|^^), (197) 



X 



where 



Ax — W X 



-2d,J 



Ay = lUr^'^^ (198) 

Defining the matrix F as G{\\ x — y \\) — ^.^^^^~1a , we will have from (196) and (197), 

\\^ y\\ 

AF(||x-y||)A* = F(A||a;-y II), (199) 

and 

A,F(|| X - y \\)Al = Fi f-^'y] ). (200) 

II a: nil y \\ 

For every arbitrary A, we now choose the points x and y such that || x \\— X~ and 

_3 

II y \\— X~ . It should be noted that in this way by varying X, we can span all the points 

1 3 

of space. From (198), we will have A-^ = A2 and A^ = A2, therefore eq. (199) turns into, 

A5F(|| x-y ||)(Ai)* = F{X\\x-y ||). (201) 
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Combining (199) and (201) and using invertibility of A, we arrive at, 

F = A5F(A*)^, (202) 
by iterating (202), we will have F = AFA*, and by rearranging, we have, 

FA* = AF. (203) 

Expanding A in terms of power of In A as A = 1 + (c? In A) J + ^^i'^^ + • • • and comparing 
both sides, we arrive at 

F{J*f^{jfF, A; = l,2,---,n-l. (204) 
Since {J^)ij = Sij+k, we will have from (204), 

Fi,j-k — Fi-k,i, (205) 

which means that on each opposite diagonal of the matrix F, all the correlations are 
equal. Moreover from FJ = JF, one obtains, 

Y,Fa5^,i+^^Y.5,,i+,Fi^, (206) 
1=1 1=1 

which means that if j = 1 and 1 < i < n, then Fj_ij = 0, or 

Fij = for j = 1 and 1 < i < n - 1. (207) 

Combining this with (205), we find that all the correlations above the opposite diagonal 
are zero. In order to find the final form of F, we use eq. (199) again, this time in 
infinitesimal form, let A = 1 + aJ + o{a^) where a = (ilnA, then from AF(a;)A* = F(Aa;), 
we have, 

dF 

d(JF + FJ')=x—. (208) 
dx 

Due to the property (205) only the last column of F should be found, therefore from (208) 
we obtain, 

dFin ,„ _ 

dF- 

x—^ = 2dFi,„_i, if i>l (209) 
dx 

which upon introducing the new variable y = 2d\iix gives, 

Fi,n = ci, F2,„ = ciy + C2, F3,„ = ^ciy^ + C2y + cs, etc. (210) 
with the recursion relations, 

dF- 

^ = Fi-ln. (211) 

dy ' ^ ^ 
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Thus we have arrived at the final form of the matrix F, which is as follows: 





( •• 








%\ 




•■ 







9i 


F = 


•• 


• 9o 


9i 


92 



\9o 



9n 



9n-l 9n I 



(212) 



where each gi is a polynomial of degree i in and Qi = . All the correlations depend 
on the n constants Ci, ■ ■ ■ , c„, which remain undetermined. We have checked (as the reader 
can check for the single n = 2 case) that inversion symmetry puts no further restrictions 
on the constants q. 

The observation that the transformation properties of the members of a Jordan cell 
are as of the formal derivative of the top field in the cell, allows one to determine the 
correlation functions of all the fields within a single or different Jordan cells, once the 
correlation function of the top fields are determined. As an example, from ordinary CFT, 
we know that conformal symmetry completely determines the three point function up to 
a constant. Let 0a, (pp and (p^ be the top fields of three Jordan cells (Aq, I), {Ap, m) and 
(A-y,n) respectively. Therefore we know that. 



'■a/37 



X 



2; II ^ '^11 y 



|A3+A^-A„' (213) 



where the constant Aajs^y in principle depends on the weights Aq, and A^. Denoting 
the second field of the cell (A^, /) by 0ai/ we will readily find from eq.(213) that. 



d 



+ 



A' 



a/37 



x-y 



|A<:«+A^-A^| 



X — Z 

dA. 



I A„+A^-A 



7-^/3 1 



y-z 



A^+A^-Ao 



a/97 



X-y 



In I 



I A^+A^-A^ I 



y- z 



X — z 



I A„+A-,-A^ I 



y- z 



A^+A^-Ao 



X-y 



X — z 



(214) 



where A^^^ = —d-^-Aap'^ is a new undetermined constant. For the correlation functions 
of fields within a single cell, one should then take the limit /3, 7 — > a in the above formula. 
It's not difficult to check that this formula satisfies all the requirements demanded by 
conformal symmetry. 



^For simplicity, we have denoted the top field by <j) and the second field by <j)i, instead of <j)i and ^2 
respectively. 
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4 -Disordered Systems & Logarithmic Conformal Field 
Theory 



• Introduction 



Consider a renormalization group transformation acting on the space of all possible 
couplings of a model, {k} [1]. The transformation has the form {k'} = R{k} where R 
depends, in general, on the specific transformation chosen, and in particular, on the length 
rescaling parameter b. Suppose there is a fixed point at {k} = {k*} the renormalization 
group equations, linearized about the fixed point, are 



k'a-k:^^T,,{k,-kl) (215) 

b 

where Tab = I k=k* ■ Denote the eigenvalues of the matrix T by A* and its left eigenvectors 



dkb 

by {0*} so that 



^</.lT„b = A>; (216) 



Now we define scahng variables Ui = J2a 4>a{ka — A;*), which are hnear combinations of the 
deviations ka — k* from the fixed point which transform multiplicatively near the fixed 
point 



Ui = j:<t>iiK-K)-j:<t>aTab{ka-k:) 

a. a,b 

= Y,)^'mb-K) = yu, (217) 

b 

It is convenient to define the quantities t/i by A* = 6^% the 7/j's are called RG eigenvalues. 
Now we can distinguish three cases, 

1- If Hi > 0, Ui is relevant, 

2- If Hi < 0, Ui is irrelevant, 

3- If Hi — 0, Ui is marginal. 



• Rule of the rescaling factor b 

Consider an infinitesimal transformation, that b = 1 + SI, with SI « 1. In this case the 
coupling transforms infinitesimally 

ka^ka+ (^)Sl + 0({Slf) (218) 
and the RG equation has the differential form 

^ = -mk}) (219) 
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where the function /?„ are called the RG beta-function. Also the matrix of derivatives 
at the fixed point is now Tab = ^ab + with eigenvalues (1 + 5iy^ ~ 1 + yi5l. 

Hence the y^'s are simply the eigenvalues of the matrix — evaluated at the zero of the 
beta-functions (because at the fixed points we have ka — > ka and beta-function must be 
zero) . 

• The perturbative RG 

When two fixed points are sufficiently close, it is then possible to deduce universal 
properties at one fixed point in terms of those at the other. Such an analysis is the basis 
of the e-expansion and many other similar techniques. It also allows us to describe the 
properties of fixed points with exactly marginal scaling variables. 

Now consider a fixed point Hamiltonian H* which is perturbed by a number of scaling 
fields, so that the partition function is [1]: 

Z = TreM-H* - Y^giY^^^^^Ur)} (220) 

i r 

Let us expand this in powers of couplings Qi, 

z = z*{l-Y9^^ <m>A 

■J Cii 



v 

ijk 

+ ■■■} (221) 

where all correlation functions are to be evaluated with respect to the fixed point hamil- 
tonian H*. We now implement the RG by changing the microscopic cut-off a ^ ba , with 
b — 1 + 51, and asking now the couplings gi should be changed in order to preserve the 
partition function Z. The length a appears in three ways in eq.(221), 

1- Explicitly , through the divisors a'^"^' 

2- Implicity , we must restrict all integrals to |rj — rj\ > a. 

3- Through the dependence on the system size L in the dimensionless ratio L/a. 
For the explicit dependence we have: 

^ (1 + Sl)''-^^gi ^gi + {d- Xi)g,5l (222) 

The implicit dependence will appear in the second order term. After changing a 
a{l + SI) we may break up the integrals, 

/ = / - / (223) 

J\ri-r2\>a{l+Sl) J\ri-r2\>a J a{l+Sl)>\ri -r2\>a 

The first term simply gives back the original contribution to Z and the second term may 
be expressed using the operator product expansion 

<l>i{ri)<pj{r2) = E Cijkiri - r2)Miri + r2)/2) (224) 
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where Cijk{ri — = |^ _j. '^^x^Xx^-xf. , (cjjfc are known as OPE coefficients) as: 



J ,, . a(l+50>|n-r-2|>a 

<4>.{{n + r,)/2)>p$§: (225) 

The integral gives a factor Sdof'Sl, where — 27r'^/^/r((i/2) ( the area of d-dimensional 
sphere). This term may then be compensated by making the change 

9k^ 9k- 1/25'd Cijk9i9jSl (226) 

Finally by changing the variable —>■ {2/ Sd)9k we find the beta-function for coupling 
as (with Uk = d — x^)'- 

^ = yk9k - Cijk9igi H (227) 

i,3 
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• Quenched Random Ferromagnets 



We consider the random bond Ising model and suppose that positions of the impurities 
are fixed , and tracing over the only magnetic degrees of freedom [5,60]. Let us describe 
these disordered systems in the continuum limit by the following Hamiltonian, 

H^H'^ + J j{r)E{r)dfir (228) 

where is the Hamiltonian of the renormalization group at fixed point describing the 
pure Ising model. The field j(r) is a quenched random variable coupled to the local energy 
density E{r-). For simplicity we assume that the random variable j(r) is a gaussian variable 
which is characterized with its two moments < j{r) >— and < j{r)j{r') >— g5{r — r'). 
The exphcit distribution of j(r) is, 

P{j{t)) = NeM-^ I J^rfd^r) (229) 

where N is normalization constant. 

We are interested in computing the average of the free energy, or quenched averaged 
correlation functions. Since the free energy is proportional to logarithm of the partition 
functions Z[j], we have to compute the average of logZ[j]. The replica method is based 
on the identity 

\ogZ\j] = hm^^^ = ^Z-[j]\n=o (230) 

n-i-o n an 

The standard procedure of averaging over disorder is to introduce replicas , i.e. , n 
identical copies of the same model 

= Tr exp{- E - / ^'^j(^) E ^a(r)}. (231) 

a=l ■' a=l 

Averaging over the disorder gives rise to the following effective replica Hamiltonian: 

n „ n 

Hn = Y:Hl-g Ea{r)E,{r)d''r. (232) 

a=l a^b 

We keep only the non-diagonal terms , since using the operator algebra of the pure system 
one can absorb the diagonal terms into H^. The replicas are now coupled via the disorder. 
The scaling dimension of coupling g is yg = d — 2x^^ and it is relevant at the pure fixed 
point \{ d > 2x%. If yg is small , it is possible to develop perturbativc RG in powers of 
these variable and we can find a random fixed point with perturbativc RG. 

To find the renormalization group equation for g we need the operator product expan- 
sion of {J2a^b EaEb) with itsclf. Siucc the replicas are decoupled in , we may evaluate 
this using the operator product expansion of with itself, the first few terms of which 
have the form 

Ea-Eb^ 5ab + cAbEa + ■■-, (233) 

where c is a coefficient whose value is fixed and universal. The coefficient c vanishes when 
d = 2. This is a consequence of the self-duality if Ising model in two dimensions [1]. Now 
we have 

{J2EaEb)-{J2EcEd)^ E EaE,E,Ed. (234) 
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The OPE of EaEb is zero (see eq.(238)) and the OPE of EcE^s is also zero, so we have 
two ways of writing the OPE. So 

X] EaEbEcEd 

= 2 ^ {5hd + c6bdE,)EaE,. (235) 

The first term is 

2 ^ 5bdEaEc 

= 2{(n -1)Y: El + (n-2)Y: E^E,}, (236) 

a=l a^b 

the second term can be written as; 

2c ^ 6bdEaEcEb 

a^byC^d 

= 2C ^ 55d(5„c + c5acEa)Eb 

a^b,c^d 

= 2C ^ dbd^acEb + 2C^ ^ SbdSacEaEb 

ayth,Cytd a^b.c^d 

= 2c(n-l)^E„ + 2c2 ^ SuSacEaEb 

= 2c(n - 1) ^ E„ + 2c' ^ (237) 

a=l a^^b 

Therefore: 

{Y.EaEb)-{Y.E,Ea) 

= (2(n-2) + 2c2)^Ea£;6 

+2{n + 2c{n - 1) ^ ^« (238) 

a=l o=l 

The renormahzation group equation for g is thus 

I = 2/;^ - (2(n - 2) + 2c')/ + 0(/, • • •)■ (239) 

Therefore if one denote Ea^b -E'a-E'b with $ so the coefficient in OPE $ • $ = 1 + 6$ H 

, is 2(n - 2) + 2c2. 

Now Consider a scahng operator (p with a scahng dimension a;^ so that 0$ = 6^0 + • • • 
and denote couphng of (p with t. Now we have 

/3, = ^j5-6/ + --- 
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= ylt-2b^tg + --- (240) 

Prom the above equations, 

/ = f (241) 

Therefore using eq.(245) we obtain: 

y. = ^lo = 4-2^'4' (242) 

and then 

d-x^ = {d-xl)-^, (243) 



so 



Now consider the case where ((> is Et = J2a and Ea = Ea — {^) J2a Ea with Ea = 0- 
The combination Et = J2a=i is a singlet ( symmetric under the permutation or the 
rephca group) and the Ea = Ea, — ^ J2b=i Eb transforms according to an (n— l)-dimensional 
representation of Sn- The fields Et and Ea have proper scaling dimensions. To find the 
scaling dimensions of new fields we should find the OPE coefficients Ef^ and Ea^, 

n n 

Et<^ = {J2Ea)iJ2E,E,) 

o=l b^c 



n 



— X] EaEbEc 
a,b^c 

= 2j2{^ab + c5abEa)E, 

a,b^c 
n n 

= 2 X SabEc + 2C X! i^ac + c5acEa)5ab 

a,bj^c a,b^c 
n 

= 2j2SbbEc + 0{c)^2{n-l)J2Ec, (245) 

b^c c 



SO — 2{n — 1). 



Ea^ 



in n 

{Ea--Y.Eb){ E EcEa) 



n 



6=1 c^d=l 



Y,EaE,Ea--{Y.E,){ E E,Ea) 



b=\ cj^d=l 



= 2j2{Sa. + cSa.Ea)Ed--{^Et>){^ E,Ed) 



n 



b=l c^d=l 
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n 1 n 

= 2Y,Ea--{2n-2){Y.E,) + --- 

di^a ^ 6=1 

1 

= -2(£;a--E^^') + --- (246) 

6=1 



SO h^^ = —2. Now we can obtain , x^^, 

2h 

XE, = x% + -^y, (247) 
we have neglected 0(c) and 0(c^ • ••) so 6 = 2{n — 2) and in the hmit of n — > 0, 



x^;, = 4 + (l+n/2)y^ (249) 
= 4„ + (1 - ^/2)y., (250) 
Now define = SA^^ and a;^^^ = 2Ag_^. 

The important observation is that the fields Et and Ea have the proper scaling dimensions 
close to n ^ as A^, = A'§l+yg/2+0{y'^) and A^ = A^^+yg/2+0(y^) respectively. It is 
clear that the singlet field Et becomes degenerate with the (n — l) operators However 
they do not form the basis of the Jordan cell for the dilatation operator. Starting from 
the following replicated Hamiltonian one can show that correlation functions calculated 
against this effective Hamiltonian will correspond to correlator averaged against the initial 
Hamiltonian with quenched disorder, 

n „ n . 

HR^Y.Ha + t d''xY,Ea{r)-g d'r ^ Ea{r)E,{r) 

a=l •' a=l ^ a^b 

< E{r) >H ^ lim < Ea{r) >repi 

< E{Q)E{r) >H ^ lim < ^„(0)^„(r) >repi 

< E(0) >H< E{r) >H ^ lim < Ea{0)E,{r) >repi a ^ b 

< E{ri)E{r2)E{r3) >h ^ hm < Ei{l)Ei{2)Ei{3) >repi 

n— >0 

< E{1)E{2) >H< E{3) >H ^ ]im < Ei{l)Ei{2)E2{3) >repi 

n—^O 

< E{1) >H< E{2) >H< E{3) >H 

^ lim < ^1(1)^2(2)^3(3) >repl (251) 
n— >0 

etc. 
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To find the logarithmic pair consider: 



< Et{0)Et{r) >= 

n(< Si(0)Ei(r) > -{n - 1) < Ei{0)E2{r) >) 

< Ea{0)Ea{r) >= 

(1 - -)(< Ei(0)E,(r) >-< E,(0)E2{r) >) 
n 

= (1 - -)S(n)r-2^^(") (252) 
n 

The above equations enable us to write the quenched averaged connected two-point 
correlation functions of energy density operator in terms of < Ei{0)Ei{r) > and < 
Ei{0)E^{r) > in the limit of n ^ as: < E{0)E{r) >^ =< Ei{0)Ei{r) > - < Ei{0)E2{r) > 
which is equal S(0)r~^^^ and it has a pure scaling behavior. However the correlation 
functions < £^i(0)£^i(r) > and < Ei{0)E2{r) > have the logarithmic singularities and 
behave as: 

< Ei{0)Ei{r) >= {A'{0) - B'{0)+B{0) 
- S(0)^ In r)r-2^^ 

<Ei(0)E2(r) >= (yl'(O) - £'(0) 

-^(0)^lnr)r-'^^ (253) 

where ^4(0) = -B(O). The prime sign in the eq. (258) means differentiating with respect 
to n. 

This means that in the limit n — > the field Ei and Ea form a basis of Jordan cell, i.e. 
their two point correlation functions behave as: < Et{0)Et{r) >= 0, < Et{0)Ea{r) >~ 
air"'^^'^ and < Ea{0)Ei,{r) >— (— 2ailnr + Da,b)r~'^^, where ai and Da,b are some con- 
stants. 

Also the ratio of quenched averaged two-point correlators of the energy density operator 
to connected one has a universal r- dependence as: 



<E(0)E(r)> <E(0)><E(r)> ,^ 
<B(0)B(r)>„ <B(0)B(r)>„ 

We note that in 2D we have dealt with two-dimensional conformal field theory, relying 
heavily on the underlying Virasoro algebra. For an extension to d dimensions one has to 
modify the representation of the Virasoro algebra to higher dimensions. We consider a 
doublet of fields (Jordan cell) 
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*=(|:). (255) 

and note that under D-dimensional conformal transformation x — > x', we have, $(x) 
$'(x') = G"'"$(x) where T is a two dimensional matrix which has Jordan form and 
G — 1 1 is the Jacobian. The fields and Ea, transforms as: 

Ei(x') = G-—Et{^) 

Eai^') = G-—{ln{G)Et{x)+Ea{x)) (256) 



This expresses that the top-field Et always transforms as an ordinary scaling operator. 
It can be verified that the correlation functions of fields E^ and Ea have the standard d- 
dimensional logarithmic conformal field theory structures. 
Let us consider the case that d — 2. Therefore we can write: 



LoEt = AeEi 
LoEa = AeEu + Ef 
LoEb — AeEj, + Et 
LqEc — AeEc + Et 

where we have used the replica symmetry. Using the above equations, it is evident that 
the dimension of difference-fields Ea — E^ with a ^ b is Ae and it transforms as an or- 
dinary operator under the scaling transformation. The important observation is that the 
individual logarithmic operator Ea do not contribute to the connected quenched aver- 
aged correlation functions. Instead the connected averaged correlation functions depend 
on the difference fields Ea — Ef, only. For instance in the following we write the con- 
nected quenched averaged 2,3 and 4-point functions of local energy density in terms of 
the difference operators explicitly, 

< E{1)E{2) >^ = ^<{Ea-E,\,){Ea-E,\,)> (257) 

< E{1)E{2)E{3) >^ = <(£;„-E,)(i) 

{Ea-Ec)i2){Ea-Eh)i3)> (258) 

< E{1)E{2)E{3)E{A) >^ = 

< {Ea — -E'b)(i)(-E'a - Ec){2){Ea - Ed)(3){Ea — -E'(,)(4) > 
~2 ^ ~ ~ Ed){2){Ec — E(i)(3){Ea — -E'6)(4) > 

< {Ea - -Eb)(i)(-Ec - Ed){2){Ea — Eb)^3){Ec - -Erf) (4) > 



40 



— ^ < {Ea — Eb)(^i^{Ea — Eb)(2){Ec — Ed)(3){Ec — Ed) (4:) 



> 



where the last equation has only 15 independent terms. This shows that quenched av- 
eraged connected correlation functions have a pure scaling behavior. Let us verify this 
results from direct calculation of quenched averaged connected 3-point correlation func- 
tion of energy density. 

We are interested in exact derivation of the various 3-point quenched averaged functions 
as < E{1)E{2)E{3) >, < E{1)E{2) >< E{3) > and < E{1) >< E{2) >< E{3) >, which 
can be written in terms of the replica correlation functions < £'i(l)£'i(2)£^i(3) >= a 
< ^1(1)^1(2)^2(3) >= b and < ^1(1)^2(2)^3(3) >= c, respectively One can derive the 
correlation functions a, b and c by means of the 3-point functions of Et and Ea as: 

< Et{l)Et{2)Et{3) >= na + 3n{n - 1)6 

n{n - 1) (n - 2)c = nAi (259) 
< Ea{l)Ea{2)Et{3) >= ma + {ni{n - 1) 



-4nl + ^{n - 1)2 + ^{n - l)(n - 2))6 



+ (--ni(n-l)(n-2) + ^(n-2)2(n-l))c 



n ' ' ' 11? 



= (1 - -)Si (260) 
n 



and finally, 



~ ~ ~ /?/ — 1 

< E„(1)E,(2)E„(3) >= [nl —)a 



Tt 13 3 
{-2>nl— -{n - l)(n - 2) + -^ni{n - l))b 



+(-^ni(n - l)(n - 2) - ^(n - l)(n - 2){n - 3))c 

(261) 

where ni = (1 — -). To derive the above equations we use the replica symmetry and 
symmetries of the various type of 3-point correlation functions under interchanging of 
positions. We note that replica symmetry leads to have < Ea{l)Et(2)Et{S) >— and 
therefore, dose not give any new relationship between a, b and c. Using the above equa- 
tions, it can be found that the correlation functions a, b and c are as following: 

SnBi - 3nCi + n^Ci + - 3Bi + 2Ci 
a = 

nBi - nCi + Ai- 3Bi + 2Ci 
^ = 2 
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Ai - Wi + 2Ci , , 

c = — 1 262 

Where Ai , Bi and Ci are pure scaling functions of variables r^j. Using the above 
equations we can show that the connected quenched averaged 3-point function behaves 
as: 

< E{l)E{2)E{'i) >c = 2c + a - 36 = Ci (263) 

which is a scaling function and confirms the observation that the logarithmic operators 
(individually) do not play any role in the connected quenched averaged correlation func- 
tions. In addition one can derive the correlation functions < £'j(l)£'j(2)i5fc(3) > for given 
i, j and k in the limit of n ^ and show that they have the following form: 

< Ei{l)Ej{2)Ek{2>) > = K-fe-AjfeA 



+ 7i.fc(4i^2-i^D]/(l,2,3) (264) 

where /(I, 2, 3) = {ri2ri2,r2z)~'^^'' , -Di = In(ri2ri3r23) and D2 = lnr23 Inris+lnris lnri2 + 
Inr23lnri2. It can also be shown that the ratio of various symmetrized 3-point functions 
to the connected one behaves asymptotically as a universal function 

^{AD2-Dl). (265) 

We generalize the above calculations to derive the various type of 4-point correlation 
functions and show that the ratio of the various disconnected to the connected one has 
the following universal asymptotic: 

~ 4 [Of - 6O2 - 3O3 - I2O4 - I8O5] (266) 
36 

where Oi = In(ri2ri3ri4r23r24r34), O2 = {In Vijlnrli + ■ ■ ■) with i ^ j k ^ I, O3 = 
(Inrjj Inrfi^ + ■ ■ ■) with i ^ j k, = (Inr^j Inr^j \nrij -!-•••) with i ^ j ^ k ^ I, and 
finally O5 = (Inrj^ Inrj^ In + ■ ■ •) with i ^ j ^ k ^ I 

One can check directly that these different types of the 3 and 4-point correlation func- 
tions have the general property of a logarithmic conformal field theory that the logarithmic 
partner can be regarded as the formal derivative of the ordinary fields (top field) with 
respect to their conformal weight. In this case, one can consider the Ea fields as the 
derivative of Et with respect to n . We emphasis that the derivative with respect to scal- 
ing weight can be written in terms of the derivative with respect to n. These properties 
enable us to calculate any N-point correlation function containing the logarithmic field 
Ea in terms of the correlation functions of the top-fields. We have shown that the indi- 
vidual logarithmic operators Ea do not have any contribution to the quenched averaged 
connected correlation functions of the energy density. We also obtain that the connected 
correlation functions can be written in terms of the difference fields which transform as 
an ordinary scaling operator. However they will play a crucial role to the disconnected 
averaged correlation functions. Also we find that the ratio of the various types of 3 and 
4-point quenched averaged correlation functions to the connected ones have a universal 
asymptotic behavior and give their explicit form. Our analysis are valid in all dimensions 
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as long as the dimension is below the upper critical dimensions. To derive the above 
results we have used the replica symmetry. Any attempt towards the breaking of this 
symmetry will change the above picture and may produces more than one logarithmic 
fields in the block and produce higher order logarithmic singularities. These results can 
be easily generalized to other problem such as polymer statistics, percolation and random 
phase sine-Gordon model etc. 
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